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Abstract 

An efficient method with no numerical diagonalization of a huge Hamiltonian matrix and 
calculation of a tedious Green’s function is proposed to acquire the exact energy spectrum 
and dynamical conductivity in a gated AA-stacking A-layer Graphene (AANLG) with the 
intrinsic spin-orbital coupling (SOC). 2N x 2N tight-binding Hamiltonian matrix, velocity 
operator and Green’s function representation of an AANLG are simultaneously reduced to N 
2x2 diagonal block matrices through a proper transformation matrix. A gated AANLG with 
intrinsic SOC is reduced to N graphene-like layers. The energy spectrum of a graphene-like 
layer is A = ± e| I . e± depends on the interlayer interaction, gated voltage and layer num¬ 
ber. £| I = + where Emg is the energy spectrum of a monolayer graphene and A is 

the magnitude of intrinsic SOC. More importantly, by inserting the diagonal block velocity 
operator and Green’s function representation in the Kubo formula, the exact dynamical con¬ 
ductivity of an AANLG is shown to be a = the sum of the dynamical conductivity 

of N graphene-like layers. The analytical form of Oj is presented and the dependence of Oj 
on e±. A, and chemical potential is clearly demonstrated. Moreover, the effect of Rashba 
SOC on the electronic properties of an AANLG is explored with the exact energy spectrum 
presented. 
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1 Introduction 


Graphene and its family members, ineluding the AA-, AB- and ABC-staeking graphenes, have 
long attraeted a lot of attention due to their striking physical properties. Graphene, a pure two- 
dimensional (2D) system, is an atomic sheet peeled off from graphite Ol |3. Carbon atoms are 
brought together and packed into a hexagonal lattice to form a graphene sheet. Such a geometrical 
structure consequently brings about a pair of low-lying linear energy bands. Electrons on the 
graphene sheet behave like the relativistic massless particles. This linear dispersion induces a 
variety of unique electronic properties, such as, electron-hole symmetry, Klein tunneling, high 
mobility at room temperature, non-zero conductivity, and anomalous quantum Hall effect (Si @1 
El 0 H [HI m Uni [m [I 2 I HSl eh ESI ES EVI. owing to the manifestation of fascinating effects, 
graphene is a promising material expected to play a vital role in technological applications, e.g., 
display screens, electric circuits, solar cells, analog electronics and photonics/optoelectronics [E8l 

Multilayer graphenes are the pile of several graphene layers held together by the van der 
Waals force. The low-energy physical properties depend strongly both on the stacking order 
and on the number of layers [|26l |2H [28l |29l |30l [3T1. The most studied multilayer graphenes 
are AB-stacking bilayer graphene [1^ |33l [34l ESI EH [371 EH EH [40l and ABC-staeking trilayer 
graphene [HTl 1^ [43l l44l . The AB-stacking bilayer graphene shows four parabolic bands around 
the Dirac points. The touch between valence and conduction bands makes the AB-stacking bi¬ 
layer graphene a zero gap semimetal. A band gap is opened by the application of a vertical 
electric field [|40l . The low energy dispersions of ABC-stacked trilayer graphene are described 
by two remarkably flat bands. The two-fold degeneracy in the band structure can be readily 
lifted by a perpendicular electrical field. Due to the progress in the fabrication and manipu¬ 
lation of graphene layers, the AA-stacking graphite and AA-stacking multilayer graphenes are 
produced ll45l 1^ . Following this, theoretical and experimental studies are conducted in order 
to explore electronic properties of AA-stacking bilayer and multilayer graphenes, e.g. infrared 
spectra, Raman spectra, Landau-level energies, absorption spectra, magneto absorption spectra, 
static polarization, and dynamical conductivity [[45l[4H[47l[48l[49l[50l[5Tl[52l[5Hl . 

The increase in the layer number, appearance of interlayer interaction, and application of 
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external field will lead to more diffieulties exploring the electronie properties of the multilayer 
graphenes. For instance, a IN x IN tight-binding (TB) Hamiltonian matrix is constructed and 
used to describe an A^-layer graphene with the nearest neighbor interactions taken into account. 
The exact diagonalization of a IN x IN Hamiltonian matrix will be cumbersome with the layer 
number N increasing. A high-rank Hamiltonian matrix gives rise to more tasks in calculation of 
the Green’s function, which is generally adopted to study the minimal and dynamical conduc¬ 
tivities. Moreover, the mirror-symmetry-breaking, caused by a vertical electric field applied to 
a multilayer graphene, also increases the difficulty in the diagonalization of Hamiltonian matrix. 
Most of studies focus on the exploration of physical properties as the layer number N <3, e. g., 
AB-stacking bilayer graphene, and ABC-stacking trilayer graphenes. Recently, the investigations 
of the dynamical conductivity of AA-stacking graphene and static polarization of AAA-stacking 
graphene have been reported [53l . A model can deal with physical properties of multilayer 
graphene in different stacking order or various layer number N under external field is inspired 
and desired. We previously presented analytical modes to exactly describe the minimal conduc¬ 
tivity of the AB-stacking multilayer graphene li54ll and exact Landau levels of the AA-stacking 
multilayer graphene lfSSlI . 

In this work, an analytical model is proposed in order to derive the dynamical conductivity 
and energy spectrum in a gated AANLG with intrinsic SOC. 2N x 2N Hamiltonian matrix of 
AANLG is decomposed into N 2x2 diagonal block matrices. An AANLG is decoupled into 
N graphene-like layers. Thus, a close form of the energy spectrum is disposed. Application of 
current analytical model to the study of the dynamical conductivity of AANLG is conducted. 
It is shown that the dynamical conductivity of an AANLG is equal to the sum of the dynami¬ 
cal conductivity of N graphene-like layers with/without intrinsic SOC. Above all, the presented 
model can efficiently and exactly give out the energy spectrum and dynamical conductivity in a 
gated AANLG with intrinsic SOC and avoid the diagonalization of a huge Hamiltonian matrix 
and calculation of associated Green’s function. 
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2 Gate-Tuned Energy Spectrum Of AANLG with Spin-Orbital 


Coupling 

Graphene is a two dimensional atomie sheet made up of earbon atoms, whieh are preeisely paeked 
in a planar hexagonal lattiee, viewed as bipartite lattiee eomposed of two interpenetrating trian¬ 
gular sublattiees. The earbon-earbon bond length is b = 1.42 and the lattiee vector is equal to 
a = \/3b. A primitive cell contains two atoms denoted as A and B. With SOC Taken into consid¬ 
eration, the Hamiltonian Hmg of a monolayer graphene is [[561 [571 


Hmg = ho + hjso + hR, 


( 1 ) 


where the first term, ho = Y,ij Oloc'lcj -|- h.c., is Hamiltonian operator of the monolayer graphene 
without SOC. cf (cj) is the creation (annihilation) operator and creates (annihilates) an electron 
at the site i (j). (Xq is the intralayer nearest-neighbor hopping between atoms A and B on the 
same graphene layer, as illustrated in Fig. 1. The second term hjso is the intrinsic spin-orbit 
interaction. The third term is the Rashba SOC, which is induced by the external perpendicular 
electric field or the interaction with substrate. The Hamiltonian operator of Rashba SOC is Hr = 
— IXr E/i,v (^^5 V X dij)zCj + h.c., where Xr is the magnitude of the Rashba SOC. S is the 
Pauli vector, the subscripts /i and v represent the spin index, and dij is the unit vector pointing 
from atom site i to its nearest neighbor j. 

Without the Rashba SOC {Xr = 0), the TB Hamiltonian matrix, spanned by periodic Bloch 
functions |A) and \B), is [[52l 


Hmg = 




ai 


ttk 


\ 

/ 


( 2 ) 


where oCk = cto/(k) = aoLj=iexp(/k ■ by), represents the three nearest neighbors on the 
same graphene plane and k is the in-plane wave vector. A is the strength of ISOC and 5^ = ±1 
represents the up or down spin, x^ — ^-l at the Dirac points K and K’. The energy dispersions are 
E = ±A/|akP + A^. 

Furthermore, the TB Hamiltonian matrix of the Hamiltonian Hmg = ho + hjso + hR, acting 
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on periodic Bloch functions |A t), |B t), i), |B |), is lISTl 


Hmg = 


A ttk 0 0 




ai 


—A —iXw 0 


0 iXn —A oCk 


0 0 


al 




O) 


where the Rashba SO interaction between |A t) and \B |) {\B t) and |A |)) is neglected because 
it is much weaker than A/? [l57ll . The analytical energy dispersions are 



^ + A^ + XjiA + . 


(4) 


2.1 Energy Spectrum Of AANLG with Intrinsic SOC 

By stacking N layer graphenes directly on each other with an interlayer distance between graphenes 
c = 3.35 A [|5^ , an AANLG is formed, as shown in Fig. 1. In the stacking direction, N atoms A 
(B) form a linear chain. The primitive unit cell contains 2N atoms, denoted as Ai, A 2 , • ■ ■Af^,Bi , B 2 , • • • , B/^,. 
The first Brillouin zone is the same as that of a graphene. The Hamiltonian of an AANLG is given 
by 


HaANLA = ^0 + Hjso + Hr , 


(5) 


where Hq is Hamiltonian operator an AANLG without SOC, and Hjso and Hr are caused by ISO 
and Rashba SO interactions. Following the discussion in the subsection above, we first close the 
Rashba effect, i. e. , Hr = 0. In the presence of an electric field, the Hamiltonian representation 
of an AANLG with intrinsic SOC, spanned by periodic Bloch functions |Ai IA 2 t): •'' t)> 

\Bi t), 1^2 t)r ■ ■ \Bn t)os a 2A X 2N matrix, reading 


H = 



y Hra 


Hab 

Hbb 


\ 


(6a) 
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where Haa,Hab,Hba, and Hbb are N x N matriees. The Hamiltonian operators Haa and Hab 
deseribe an A^-site linear ehain with the intrinsie SOC subjeeted to a parallel eleetrie field. 

The reeent researches [[5^ l60l , using ab intial calculation and TB method, show that in the 
Bernal stacking bilayer graphene and ABC stacking trilayer graphene, many intralayer and in¬ 
terlayer SOC parameters are included to fully describe the spin-orbital interactions in the TB 
model. The numerical fitting of TB parameters to the results of ab initial calculation exhibits 
that the magnitudes of intrinsic (A) and Rashba spin-orbital interactions (Xr) are layer-position- 
dependent. What’s more, the strength of the interlayer intrinsic SOC is much weaker than that 
of the intralayer intrinsic SOC. Accordingly, we neglect the interlayer intrinsic spin-orbital inter¬ 
actions in our case and then take the intralayer intrinsic and Rashba spin-orbital interactions into 
account. Thus, the two 4x4 Haa and Hrb Hamiltonian matrices of the AA-stacking quad-layer 
graphene, for instance, are expressed as the following. 


Haa = 


3y/2-fAi cti 
cti V /2 + A 2 


0 

0 


OCl 

0 


0 

OCi 

-V/ 2 +A 3 

OCi 


0 ^ 
0 
Cti 

-3y/2 + A4 


m 


Hbb = 


3y/2-Ai 

OCl 

0 

0 


OCl 

V/ 2 -A 2 

OCl 

0 


0 

OCl 

-V/ 2 -A 3 


0 ^ 
0 
cti 

-3y/2-A4 y 


(6c) 


where Aj (j = 1,2,3 and 4) is the intralayer intrinsic spin-orbital interaction of the y-layer 
graphene. V = \e\Fc is the effect electric potential difference between the adjacent layers caused 
by the external electric field. Then, V is denoted as the gate voltage. In multilayer graphenes, the 
potential drop between two adjacent layers might be affected by the screening processes lIMll . For 
simplicity, we also assume that the potential drop is the same on each graphene layer. The inter¬ 
layer hopping parameter, ai, couples the two A (or B) atoms from the two adjacent layers [Fig. 
1]. as, the interlayer interaction between atoms A and B from the two adjacent layers, results 
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in a weak electron-hole asymmetry in an AANLG f[55l l58l . The values of the hopping integrals 
are ai = 0.361 eV and 053 = —0.032 eV ll5^ . Only the main interlayer interaction ai is taken 
into consideration because of a^, a\. The matrix element Hab, resulting from the intralayer 

interaction, reads 


hab= m, = 


^ ak 0 0 0 ^ 

0 ock 0 0 

0 0 ock 0 

0 0 0 ttk y 


( 6 d) 


The layer-position-dependent intralayer intrinsic spin-orbital interaction and Rashba SOC 
destroy the inversion symmetry of AANGL. The breaking of the inversion symmetry complicates 
the analysis and discussion. For simplicity, the intralayer intrinsic spin-orbital interaction are 
assumed to be independent of the vertical positions; that is, Aj = A and Xrj = Xr. Notably, the 
Haa (Hbb) is the sum of two matrices Hy and Hjso and reads 


Haa = Hv + Hjso = Hv + A1, (7a) 

Hbb = Hv+Hiso = Hv-A1, (7b) 

where 1 is a A x A identity matrix. Hy describes an A-site linear chain without SOC subjected to 
a parallel electric field. Haa (Hbb) is commute with Hy -, that is, Haa (Hbb) and Hy share the same 
eigenfunctions. The eigenenergy £j and associated eigenfunction \Sj) of Hyaxe easily obtained 
through the diagonalization of the eigenvalue equation [l55ll 

Hy\Sj)=£j\Sj), (8) 


where 7 = 1,2, • • • , A. The transpose of \Sj) is |Sy)^ = Sj^ 2 , Sj, 3 , ■ ■ ■ Sj,N) and the component 

Sjj, is the site amplitude of atom A or B located at the /th layer. 

With column vectors IS 7 ), the eigenfunctions, an A x A unitary transformation matrix Uy = 
(|Si), IS 2 ), ■ ■ • JSa?)) is then constructed and used to diagonalize Hy, i. e., UyHyUy = £jl. 
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where 1 is a unit matrix. The eigen-energies of Haa and Hbb are £j + A and £j — A after the 
diagonalization of Eqs. ( fTa] ) and ( |7b| ), respeetively . 

To aequire the energy speetrum of AANLG, a 2N x 2N unitary transformation matrix 


U 


( - \ 
‘ Uy Q ^ 


V 


0 Uy 


(9) 


/ 


is built to transform the Hamiltonian matrix (Eq. (6a) ) into a simple form. After the unitary 
transformation, a redueed matrix Jifred = U^//U has the form 


^^red 


(e; + A)l ockl 
^ a*l {£j-A)l ^ 


( 10 ) 


where 1 is A x A unit matrix. Then, the reduced Hamiltonian matrix can be rearranged into 
block diagonal form, J^red — © =9^ © • • • © where each 2x2 block diagonal matrix M’j 

is expressed as follows 


/ 


— 


Cj + A ctk 




al 


£j-A 


( 11 ) 


That is to say, an AANEG can be decomposed into A subsystems, The exact energy spectrum 

of each subsystem is 

= Cj© + = e^±e||, 

where £; (= £±) depends on the magnitude of interlayer interaction, gated voltage and layer num¬ 
ber. £|| = y/lctkP + A^ is the energy spectrum of a monolayer graphene with SOC. 

Around the Dirac point A, the diagonal block is for k = K + q 


/ 


_ 

— 


£/ + A 




^ -hvpiqx-iqy) 


£j - A 




/ 


( 12 ) 


and hyp = jOob is the Eermi velocity. The low-lying energy dispersions associated with are 
= £j ± ^/\hvFq\^ + A'^, where \q\ = ^ql + qj. 


8 
















2.2 Energy Spectrum Of AANLG with Intrinsic and Rashba spin-orbital 


Interactions 


If we take both the intrinsie and Rashba spin-orbital interaetions into eonsideration, TB Hamil¬ 
tonian matrix of an AANLG subjeet to a perpendieular eleetrie field, aeting on periodie Bloeh 
funetions |Ai t),|Ai ;),|Bi i), ■■■, |A^-1), 1), l^j I), l^,/i), ■■■, l^yv t), t), I 

), \Bm I), is a 4A X 4N Hermitian matrix and expressed as follows 


^ Hi Ht 0 . 0 ^ 

Ht H2 Ht ■■ : 0 


HAANLG = 


0 Ht H^ Ht 0 
0 0 


0 


V 


0 


■ ■. ■ •. ■ •. Ht 

0 ■■■ Ht Hn j 


(13) 


where Hj and Ht are 4x4 bloeks. The off-diagonal bloek Ht = ail originates in the main 
interlayer interaetion ai. The diagonal bloek Hj = Vjl +Hmg is the Hamiltonian matrix of the 
y-layer graphene in the presenee of the gated potential Vj, whieh has the form 


^ Vj + A ttk 0 0 ^ 

a* Vj-A -Hr 0 

0 IXr Vj — A aji 

^ 0 0 al Vj-V A j 


(14) 


It is easy to diagonalize the bloek Hj = Vj\ 4 -Hmg through a 4 x 4 unitary transformation matrix 

U, whieh transforms Hmg into a diagonal matrix, .i. e., = Jzag(A++, A-|__, A_, A_+) 

(Eq. (Q). The eigenvalues of Hj are Ay = Vy -|-A±±. Then, a AN x AN unitary transformation 
matrix U = diagiV, U,..., U) is eonstrueted and used to transform Haanlg into a diagonal bloek 

form. After the operation, we obtain Haanlg^ = H++ ©//+- ®H _©//_+, where ( 

T] = ±, = ±) is an A X A matrix. 

We take the AA-staeking trilayer graphene as a study model. The Hamiltonian matrix Haatlg 
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and unitary transform matrix U are 


^ Hi Ht 0 ^ 


Haatlg = 


V 


Ht H2 Ht 
0 Ht H 3 


^ V 0 0 ^ 


and U = 


V 


0 U 0 

0 0 u 


(15) 


After the operation Haatlg^, Haatlg is arranged into the block diagonal form Haatlg = 
Ht+ ©//+- ®H _©//_+, where each is 3 x 3 matrix and -\-Hy. The latter 


term Hy, Eq. (7a), describes an A(= 3)-site linear chain without SOC subjected to a parallel 


electric field. The energy spectrum related to are E = A^^ + = e|| + e^, where A^^ = e|| 


and =->/y2 + 2ai^ 0, or 


3 Electronic Properties and Discussions 

The energy dispersions of an AANLG with SOC in the presence of the gated potential are easily 
obtained through the calculation of energy spectrum of each subsystem by using analytical for¬ 
mula Ej = + e||. For example, the energy spectrum of the A A-stacking bilayer graphene are 

described as E = al + v^/A + where A±± = ± ^IctkP + Since 

the energy dispersions are symmetry about E = 0, only the energy spectrum E > 0 are shown 
in Fig. 2. In the absence of the gates potential and SOC, the energy dispersions around the 
Dirac point K illustrate one pair of linear bands crossing at F = Ofi (dashed curves in the in¬ 
set). The gated potential V = OAai shifts the linear bands upward (red solid curves in the inset). 
In the AB-stacking bilayer graphene, the intrinsic SOC parameter is A ~ 10^^ eV and hence 
A/cti ~ 10~"^ ll59l . For convenience of numerical analysis, we use A/cti =0.1 in this work. The 
analytical model and numerical results are relevant and applicable to the exploration of physical 
properties in multilayer graphene-like systems. The inclusion of the intrinsic SOC A = O.ltti 
changes the linear bands (red solid curves) into the parabolic bands (blue solid curves), which 
are described by F = ij -|-v2/4 + ^/\hvFq\^ + The maximum (minimum) of the parabolic 
band, located at the Dirac point F, is F = (a/1.16 — 0.1)ai (F = (a/1.16 + 0.1)ai). The green 
curves are the energy spectrum of the gated AA-stacking bilayer graphene with the Rashba SO 
interaction Xr = 0.05 tti. The Rashba SO interaction Xr destroys the degeneracy of the linear 
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bands and produces four parabolic bands (green curves). The middle two parabolic bands touch 
each other at £ = Vl.lbai. 

The energy dispersions of the AA-stacking trilayer graphene are also evaluated with formula 
E = A±± + ej, where Ci = + la\, £2 = 0, and £3 — — \Jv^ + 2al. Without SOC, the energy 

spectrum E 2 = A±± + £2 are independent of the magnitude of gated potential, as shown by the 
black dashed and red solid curves in the inset of Fig. 3. Two linear bands cross over atE = 0. The 
intrinsic SOC A changes the linear into the parabolic bands, as illustrated by dashed blue curves, 
which is simulated by E 2 = + A^, The maximum (minimum) of the parabolic band is 

determined by the strength of A = O.lcti. There are four parabolic bands after the inclusion of 
the Rashba SOC (green cures). The middle two parabolic bands do not touch at £ = 0 due to the 
intrinsic SOC. 


4 Green’s function and Velocity operator 


After the introduction of the Rashba effect, as illustrated in the section above, the AN x AN Hamil¬ 
tonian matrix can be divided into four N xN diagonal blocks. This would complicate the discus¬ 
sion and block us to pursue a simple analytical form of the conductivity of an AANLG. We, then, 
switch off the Rashba effect and consider the intrinsic SOC (ISOC) alone in the following work. 
Now, we exhibit that the Green’s function and velocity operator associated with an AANLG can 
be transformed into the diagonal block matrices. With the Hamiltonian matrix H, it is straight¬ 
forward to calculate the Green’s function through G{z) = The larger Hamiltonian matrix 
gives rise to more complex tasks in calculation of the inverse matrix of zl — H. To reduce the task, 
we use the unitary operator U, which causes J^red = U^//U, to transform the Green’s function. 
After the operation, we have ^ = U^GU = ^red is ^ block diagonal matrix 

and so does {zl — J^red)- Now, the Green’s function ^ is also in a block diagonal form; that is, 
^ © ^2 © • • • © ^N- Moreover, each sub-Green’s function is and it is a 2 x 2 

matrix. 



/ 

V 


z-£j-A 


-az 


-Ofk 


\ 


z-£j+A 


/ 


(16) 
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with the corresponding elements 


/ \ _ z—fij+A 

- (z-£^)2^(A2+|akP) 

^hl 2 (z) = (z-e^.)^-(A2+|akP) 
^ 2,21 (z) = (^-e^) 2 _'(A 2 +|akP) 


/ \ _ ^ ^ 

^2,22 iz; - (2_e.)2_(A2+|akP) 


1 + 


A/A 2 +|akp. 


z-£ 

a.! 


7?+ 


2 v/A 2 + |akP 

«k 

2 v/A 2 + |ak |2 


+ ■ 

1 


-y/A^+lakP- 


Z-Ej,+ 

1 






z-£;,+ z- 




A/A 2 +|akP- 


z-£, 


+ 


1 + 


A/A^ + IttkP- 


■Ej,- 


(17) 


The velocity operator, V = is approximated as the derivative of the Hamiltonian with 


respect to the momentum hk, based on the gradient approximation. According to the Eq. ( 6 a), 
the velocity matrix related to AANLG is 


/ 


dtlAA dHAE 
<9k 5k 




dHsA dHsB 

\ 5k 5k / 


(18) 


Since Haa and Hbb are independent of the wave vector k, and are equal to zero. As a 
result, the velocity matrix is 






Vfl 



(19) 


here 1 is an A x A identical matrix and Vf = is the Fermi velocity. After the action of the 
transformation matrix, Y = U^VxU, the transferred velocity matrix is a diagonal block matrix, 
which is in the form of 7^ = © ^2 © • • • © %!■ Each 7^-, a 2 x 2 matrix, is 


y,= 


f 0 \ 

0 Vf 

\ V 


( 20 ) 


The unitary transformation matrix U, diagonal block Green’s function representation and 
velocity operator Y are now utilized to derive the analytical form of the dynamical conductivity 
of AANEG. 
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5 Dynamical Conductivity of AANLG 


The finite frequeney eonduetivity is studied by using the Kubo formula. The eonduetivity is 
written in terms of the imaginary part of retarded eurrent-eurrent eorrelation funetion divided 
by frequeney Q. as Oap = ^ ^ where is also referred to as the polarization 

funetion. Furthermore, the polarization funetion ean then be written in the bubble approximation 
as 




l(On 


cfk 

{iTzf 


Tr[Vo; ^p 


( 21 ) 


where Va is the veloeity operator in the direetion a =xor:y and G{iCOm k) is the Green’s funetion. 
With the speetral funetion representation 


Gnmi.Z) 


dcoAnJco^) 


'-oc 2n z — CO 


/ ’ 


( 22 ) 


the real part of the eonduetivity, at the zero temperature T = 0, is expressed as 


, , f°° doo . , . , 

^ap (^) = ^ ^ -^)-f(co + Q- At)] 

/■ d^k 


X 


(iTT) 


,Tr[VaA {(0 + kl,k)VpA{(0,k)], 


(23) 


where f{x) = l/[exp(x/r) + 1] is the Fermi funetion and /i is the ehemieal potential. Follow¬ 
ing the aforementioned method, the AC eonduetivity for AANLG ean be direetly ealeulated by 
putting the 2N x 2N Green’s funetion representation (or speetral funetion representation) and 


2N X 2N veloeity operator in Eq. ( [23] ). The larger the Green’s funetion (or speetral funetion 
representation) is, the more ealeulation tasks are. 

To make less eomplex, we first utilize the relation, TrM = Tr-[U^MU], invariant of the traee of 
a matrix (or operator) under a unitary transformation. Then, with a proper unitary transformation 
matrix, both the Green’s funetion (or speetral funetion representation) and veloeity operator Va 
are redueed to the diagonal bloek matriees. As a result, the analytieal form of the real part of 
the eonduetivity of AANLG ean be easily aeeessible. The details are as below. First, by setting 
M = Va G((0-ff2,k) Vp G{co,k), the traee of M is TrM = Tr[C/ly« G Vp GU]. Then, insert the 
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identical matrix UU^ = I between the velocity operator Vjg and Green’s function G, and the result 


TrM = Tr[ra ^ Tr^ 


is acquired. = U^GU is the unitary transformation of G and it is related to the spectral function 
representation in such a manner: 



d(0 S2/,nn{(0') 
2 n z-(o' 


(24) 


Ya = U^VaU is the unitary transformation of Va- Thirdly, after the unitary transformation, both 

the y and £/ are the diagonal block matrices. That is to say, the operator ^ 

also a block diagonal matrix, such as,^ = © -^3 © • ■ ■ • Each — yj^a yj,p is 

a two by two matrix. Finally, the relation Tr{^) = Tr{^i) + Tr{^ 2 ) + Tr{^ 2 ) H-is used 

to obtain the AC conductivity for AANLG 


^ap(^) = 2^J + X j ^2^)2 [Tr-^1 + Tr./^2 + Tr./^3 H-], 

= (25b) 

7=1 


It is shown that the AC conductivity of an AANLG is equal to the summation of the AC conduc¬ 
tivity of each subsystem and Oj^^p (^) of o^oh graphene-like layer is 



(25c) 


The AC conductivity of the AANLG can be analytically specified. As the 2N x 2N Hamil¬ 
tonian is decomposed into N 2x2 reduced Hamiltonian matrices, the effective Hamiltonian of 
each subsystem is described as Eq. (13). Furthermore, the Green’s function, spectral function 
representation and velocity operator associated with each subsystem are 2x2 matrices. Thus, 
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the analytical form of AC conductivity of each graphene-like layer without ISOC is [15^ 

~ ^intra Winter 1 without ISOC, (26a) 

Ointra = 4c7o5(0) |jU - £;|0( |jU - £; |), 

Ointer = Oq 0(0 - 2|£y - ^|), (26b) 


where intra and inter represent the contributions resulting from the intraband and interband tran¬ 
sitions, respectively. With the ISOC taken into consideration, the analytical form of AC conduc¬ 
tivity of each subsystem reads 


o-;,„(0) = cj(0,|£;-jU|), with ISOC, 


(27a) 


where < 7 ( 0 , \ej —/i|) is the conductivity for massive Dirac particles [l62ll6^ . and it is expressed 


as 


=4^s(a)0(r-A) + 


/2A\2i 

l + (^) j0[(O-2max(T,A))]. 


(27b) 


<70 T 

The dependence of Oj^xxi'^) on the chemical potential /i, Sj (or £^), and strength of ISOC is 
clearly revealed through the afore-presented formula. 

The numerically calculated conductivity (7;fx(^) and the associated conductivity of each sub¬ 
system ay^xx(O) (denoted as sub-conductivity) of the AA-stacking trilayer graphene (TLG) are 
presented in Fig. 4. Both the intraband and interband transitions contribute to AC conductivity. A 
delta peak at frequency 0 = 0 , caused by the intraband transition, is not shown here; that is, only 
the conductivity resulting from the interband transitions is shown. The conductivity (7i xx(^2), 
< 72 ,xx(^) and <^ 3 ,xxi^) of tho subsystems are illustrated in the dashed curves. AC conductivity 
of TLG are presented in the solid curve, which is < 7 x^( 0 ) = Yj=i tho superposition of 

the AC conductivity of subsystems. According to Eqs. (26b) and ( 27b| ), the profile of each sub¬ 
conductivity Oj^xxi^) = <7o — 2|£y — /i|) is governed by the step function 0(|£; — /i|). £j 


related to TLG are £i = — y 2ci;2 -|- £2 = 0, or £3 = ylaf + V'^. In the absence of the gated 

potential (V = 0) and at /i = 0, both Oi^xx and Os^xx show the absorption edge at = a/2. 
< 72 ,xx(f 2 ) contributes a constant background conductivity, which is equal to Oq (dashed curves 
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in Fig. 4(a)). AC conductivity of TLG (solid curve Fig. 4(a)) at high frequency is equal to a 
constant value, three times of do • At /i = 0.1 ai, (Ji^xx, <y 2 ,xx and as^xx in Fig. 3(b)) show step 
edges at frequencies Q/ai = \/l + 0.1, 0.2 and x/l + 0.1. As a result, there are three steps in the 
AC conductivity (solid curve in Fig. 4(b)). In the application of the gated potential V = 0.4ai, 
the absorption edges of ai^xx and as^xx occur at Q = (V2.16 + 0.1)ai and Q = (^2.16 —0.1)ai 
(Fig. 4 (c)). The intrinsic SOC (A = O.lcti) enhances the strength of sub-conductivity <72,xx in the 
region 0.2 < Q/ai < 0.6. For comparison, the solid curves in Figs. 4(a)-4(d) are plotted in Fig. 
4(e). The characteristics of <7xx(^) of TLG are dependent on /i, V, and strength of ISOC. 

The alternation of layer number N has a great influence on AC conductivity of an AANLG. 
Figure 5 displays <7xx(^) and the associated Gj,xxi^) of tho AA-stacking quad -layer graphene 
(QLG) . At high frequency, AC conductivity of QLG illustrates a constant value, which is equal 
to four times of Co- There are two steps in the AC conductivity of QLG at V = 0 (the solid blue 
curve in Fig. 5(a)). The location of the absorption edge of each sub-conductivity Oj,xxi^) is 
controlled by the step function 0(|ej —/i|). The gated-potential-dependent energy dispersions Bj 
related to QLG are e±± = ±\/b± — C, where B = ^ ^ and C = af + 

In the absence of the gated potential (V = 0), e++ = —£-+ = 2.618ai and £+_ = —e_= 

0.382ai. The first and second absorption edges appear at Q./ai = 0.76 and Q/ai = 5.2 (Fig. 
5(a)). Oi,xxi^) (<73,«(f^)) is identical to 02 ,xxi^) ((y 4 ,xxi^))- At /i = 0.3ai, absorption edges 
occur at Q/cti ^ 0.2, 1.4, 4.6, and 5.8. AC conductivity features four steps (the solid cyan curve 
in Fig. 5(b)). The gated potential V = 0.3ai modifies £++ = 2.881ai and £+- = 0.344ai and 
changes the locations of step edges (Fig. 5(c)). The intrinsic SOC A = O.lcti enhances the weight 
of AC conductivity around Q/cti 0.2. The solid curves in Figs. 5(a)-5(d) are plotted in Fig. 
5(e) to illustrate that the effects caused by the alteration of N, ji, V, and strength of ISOC on 
Oxxi^) of an AANLG are easily and clearly revealed through the analytical formula. 

6 Conclusions 

In this work, we propose an analytical model to derive the exact energy spectrum and dynamical 
conductivity in an AANLG in the presence of a bias voltage and spin-orbital coupling at the 
same footing. First, a proper transformation matrix is built and used to transform the 2N x 2N 
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tight-binding Hamiltonian matrix of an AANLG into N 2x2 diagonal block matrices. Then, an 
AANLG is redueed to N graphene-like layers. Thus, the exaet energy speetrum of a graphene¬ 
like layer is £ = ± e||. e^, the effeetive on-site energy of the graphene-like layer, is eontrolled 

by the interlayer interaetion, gated potential, and layer number. e|| is the energy speetrum of a 
monolayer graphene with SOC. Furthermore, we analytieally study the dynamieal eonduetivity 
of an AANLG, whieh is shown to be the sum of the dynamieal eonduetivity of N graphene-like 
layers with/without SOC. The dependenee of the dynamieal eonduetivity of eaeh graphene-like 
layer on the ehemieal potential, e^, and the strength of SOC is elearly demonstrated. Above all, 
our model ean effieiently and exaetly aequire the energy speetrum and dynamieal eonduetivity in 
a gated AANLG with SOC. 
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Figure Captions 


FIG. 1. The geometric structure of the AA-stacking multilayer graphene and the intralayer and 
interlayer interactions. 

FIG. 2. Calculated energy dispersions around the Dirac point K of the AA-stacking bilayer graphene 
for different gateed potential V, intrinsic SOC A, and Rashba SOC Xr. The dashed curves in 
the inset: (V, A, = (0,0,0); red solid curves: (0.4,0,0)ai; blue solid curves: (0.4,0.1,0)ai 
; green solid curves: (0.4,0.0,0.05) cti. 

FIG. 3. Calculated energy dispersions of the AA-stacking trilayer graphene for different V, A, and 
Xr. The dashed curves in the inset: (y,A,Ai;) = (0,0,0); red solid curves: (0.4,0,0)ai; 
blue solid curves: (0.4,0.1,0)ci;i; green solid curves: (0.4,0.1,0.05)ci;i. 

FIG. 4. AC conductivity Oxx (solid curves) and sub-conductivity (Ci, 02 , < 73 ) (dashed curves) of 
the AA-stacking trilayer graphene for different V, A, and /i are presented, (a) (y,A,/i) = 

(0,0,0). (b) (y,A,At) = (0,0,0.1)ai. (c) (V^A.ii) = (0.4,0,0.1)ai. (d) (y,A,At) = (0.4,0.1,0.1)ai. 
The solid curves in (a)-(d) are plotted in (e). 

FIG. 5. Same plot as Fig. 3 but for AA-stacking quad-layer graphene, (a) (y, A,/i) = (0,0,0). (b) 

(y,A,At) = (0,0,0.3)ai. (c) (y,A,^) = (0.3,0,0.3)ai. (d) (V^A.ii) = (0.3,0.1,0.3)ai. 

The solid curves in (a)-(d) are displayed in (e). 
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